Supersymmetry and its spontaneous breaking in the random field Ising model 



Matthieu Tissier and Gilles Tarjus 
LPTMC, CNRS-UMR 7600, Universite Pierre et Mane Curie, 
boite 121, 4 PI- Jussieu, 75252 Paris cedex 05, France 



o 

(N 



•4— > 

in 

+-> 



-a 
c 

o 

> 
(N 

oo 

o 



X 



We provide a resolution of one of the long-standing puzzles in the theory of disordered systems. 
By reformulating the functional renormalization group (FRG) for the critical behavior of the random 
field Ising model in a superfield formalism, we are able to follow the associated supersymmetry and 
its spontaneous breaking along the FRG flow. Breaking is shown to occur below a critical dimension 
don — 5.1 and leads to a breakdown of the "dimensional reduction" property. We compute the 
critical exponents as a function of dimension and give evidence that scaling is described by three 
independent exponents. 
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The random field Ising model (RFIM) is one of the 
archetypal disordered systems^ and, ever since the sem- 
inal work of Imry and Map, its long-distance proper- 
ties have provided puzzles that largely remain pending. 
About thirty years ago, Parisi and Sourlas in a beau- 
tiful 2-page letter— related the critical behavior of an 
Ising ferromagnet coupled to a random magnetic field 
to a supersymmetric scalar field theory and showed that 
the supersymmetry (SUSY) leads to a "dimensional re- 
duction" (DR) property by which the behavior of the 
RFIM in d dimensions is identical to that of the sys- 
tem without disorder in d — 2. It is well established that 
the SUSY construction for the RFIM and the associated 
DR property actually break down in low dimensions j2i2iii 
but the issue has not been satisfactorily settled. In a 
recent work, we showed by means of a nonperturbative 
functional renormalization group (NP-FRG) approach 6,7 
that the breakdown of DR is related to the appearance 
of a nonanalytic dependence of the effective action in 
the dimcnsionless fields. Physically, this arises from the 
presence of rare collective events known as "avalanches" 
or "shocks". However, our formalism could not at all 
address the question of SUSY and its breaking. The pur- 
pose of the present letter is to fill this gap and provide a 
complete picture of the critical behavior of the RFIM. 

Our starting point is the RFIM field-theoretical de- 
scription in terms of a scalar field 4>{x) in a rf-dimcnsional 
space and a bare action S[<p;h] given by 
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where f x = J d d x, U B {4>) = (r/2)^ 2 + (w/4!)0 4 , and 
h(x) is a random magnetic field that is taken from a 
Gaussian distribution with a zero mean and a variance 
h(x)h(y) — Ab 5^(x — y). Taking advantage of the fact 
that, at long-distance, the thermal fluctuations are neg- 
ligible compared to those induced by disorder (formally, 
the critical behavior is controlled by a zero-temperature 
fixed pointi) , one can focus on the solution of the stochas- 
tic field equation^ 
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Provided the solution is unique, the correlation func- 
tions of the cj> field are then obtained from appropri- 
ate derivatives of a generating functional that can be 
built through standard field-theoretical techniques.— One 
first introduces auxiliary fields, a bosonic "response" 
field (f>(x) and two fermionic "ghost" fields tp(x) and 
ip(x), as well as linearly coupled sources, and one ex- 
plicitly performs the average over the Gaussian disor- 
der. After constructing a superspace by adding to the d- 
dimensional Euclidean space with coordinates x — {x^} 
two anti-commuting Grassmann coordinates 6, (satisfy- 
ing O 2 = O 2 = 00 + 08 = 0), the resulting functional can 
be cast in a superfield formalism^ 



Z[J] = J P$exp + J J(x)$(x] 
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where we have introduced the superfield <&(x) = 4>{x) + 
9ip(x) + \p(x)8 + 08cj)(x), the supersource J{x) — J(x) + 
0K(x) + K(x)0 + QOJ(x), and the superlaplacian A ss = 
d 2 , + ABdgdg; x — (x, 0, 0) denotes the coordinates in 
superspace and f = j d d xdOdO. The </>-field correlation 
functions are obtained by functional derivatives of Z[J\ 
with respect to J (evaluated for K = K = J = 0). 

The action S ss is invariant under a large group of 
transformations that mix bosonic and fermionic fields 
(hence the name, SUSY): translations and rotations in 
the d-dimensional Euclidean and 2-dimensional Grass- 
mannian subspaces and "superrotations" that preserve 
the superdistance, x 2 = x 2 + -^00. These superrotations 
can be represented by the generators = x^dg+^-Od^ 
and Q M = -x^dg + ^0<V The presence of this SUSY, 
more precisely of the superrotations, was shown to lead to 
DRii2 One knows, however, that the whole formal con- 
struction collapses when the stochastic field equation has 
more than one solution, which is the case in the region 
of interest £ 
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We propose a resolution of the above problem that al- 
lows one to study the SUSY and its spontaneous break- 
ing. To this end we upgrade our NP-FRG approach^ 
to a superfield formulation. The key points involve (i) 
adding an infrared (IR) regulator that enforces a progres- 
sive account of the fluctuations of both the <f> field and 
the disorder while ensuring that the initial condition of 
the RG flow satisfies the SUSY, (ii) considering copies of 
the original disordered system, which gives access to the 
full functional field dependence of the renormalized cu- 
mulants of the disorder, (iii) enforcing that a single solu- 
tion of the stochastic field equation (the ground state) is 
taken into account for each copy, and (iv) using the Ward- 
Takahashi (WT) identities associated with the SUSY to 
ensure that neither the regulator nor the approximations 
explicitly break the SUSY. We stress that the introduc- 
tion of copies (or replicas) is necessary to describe nonan- 
alyticities in the renormalized cumulants stemming from 
the occurence of " avalanches" , even in a superfield for- 
malism. This, unusual, combined use of supersymmetric 
formalism and replicas is central for overcoming the flaws 
of the Parisi-Sourlas construction. 

Extending our previous work to the superfield theory, 
we introduce a generating functional of the correlation 
functions at the running scale k for an arbitrary number 
n of copies of the system (coupled to the same random 
field but submitted to different external sources), 

Z k [{J a }] = [f[ 2?$ a exp [ - AS k [{<S> a }} 
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The n-copy action in the above equation is invariant un- 
der the S n permutational symmetry and a global Z 2 sym- 
metry as well as the translations and rotations in the Eu- 
clidean and Grassmannian subspaces. These symmetries 
are then shared by a quadratic regulator of the form 

ASk = \ Yu ii ^a(x)n k , ab {x,x')^ h (x'), (6) 
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with 



TZk,ab(x,x') = SabSge'Rk(\x - x'\) + R k (\x - x'\), (7) 

where 5gg, = (8 - 8')(6 - 8'); R k and R k are IR cutoff 
functions. The regulator is chosen such that it suppresses 
the integration over modes with momentum |g| <C 
and both functions R k and R k go to zero when k — > 0. In 
addition to the above mentioned symmetries, ^-[{.Ta}] in 
Eq. (5) is invariant under the superrotations when the 
sources J a , K a , K a are set to zero for all copies but one. 
The theory then reduces to a 1-copy problem. (AS k can 



be made explicitly invariant under the same conditions 
by choosing lZ kyaa to be a function of the superlapla- 

cian Ass only; as a result, R k (q 2 ) = -A B d q 2R k (q 2 ), 
where q denotes the momentum in Euclidean space.) The 
regularization ensures that the modified stochastic field 
equation has a unique solution at the microscopic scale A 
and guarantees that the theory is indeed supersymmetric 
when k = A. 

The central quantity of our NP-FRG approach is the 
effective average action,°A2 which is the generating func- 
tional of the "1PI vertices"— and is obtained from log Z k 
by a (modified) Legendre transform, 

Tfc[{$a}] = -log.Zfc[{.7„}]+5^ / J a (x)^ a {x)-AS k [{^ a }} 

(8) 

with <J> a (x) = 5 (log Z k )/ 5 J a (x). Its flow with the IR 
scale k is described by an exact RG equation (ERGE), 10 

d t r k {{<S> a }} = ^T r {d t TZ k V k [{<S> a }}} , (9) 

where t = log(fc/A) and the trace involves summing over 
copy indices and integrating over superspace; the mod- 
ified propagator "Pfe.afcfei)^) ^ s the (operator) inverse 



of r[ 2 ^ + lZ k where T^ 1 {{$ a }] is the second functional 
derivative of the effective average action with respect to 
the superfields & a (x). (Here, superscripts in parentheses 
always denote functional derivatives with respect to the 
field arguments.) 

If for each copy a single solution of the stochastic field 
equation is taken into account, it is easily shown that 
both the random generating functional has a property 
of "Grassmannian ultralocality" , which translates in the 
fact that both log Z k [{J a }] and the effective average ac- 
tion T k [{$a}] have a linked expansion in sums over copies 
and integrals over Grassmann coordinates: 
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where & a (8) denotes a superfield at the Grassmann coor- 
dinates 9, 9; (other coordinates are omitted for simplic- 
ity) and only depends on p superfields at p Grass- 
mannian "locations" (on the other hand the dependence 
on Euclidean coordinates is completely general). The 
functional T k . p [(pi, 4> p ] is related to the pth cumu- 
lant of the renormalized disorder and, correspondingly, 
rjj. 1 ^" 1 ^]., 4>p] is related to the pth cumulant of the 
renormalized random fields By inserting the above for- 
mula in Eq. (9), taking derivatives with respect to the 
superfields, and restricting the superfields to configura- 
tions <& a (x) = (j) a (x), one obtains a hierarchy of cou- 
pled ERGE's for the "cumulants" r^'^^i, <f> p ]. It 
is worth stressing that to obtain the flow equation for 



T^"p l \<f>i, •••! 4>p] with its full functional dependence on 
the p field arguments, one needs to consider at least p 
copies. Formally, the whole hierarchy of flow equations 
for the cumulants can thus be obtained by considering 
an arbitrary large number of copies. As an illustration, 
the ERGE for the first cumulant reads 



d t T ka [&] = ^J^d t R k (g 2 )P k]q - q [0x] 

+ d t R k {q 2 )P k - q - q [(j)i,<t>i\ 



(11) 



where P k [<j>] 



lfj[</>]+#fc) and Pk[<hM 



Pkifa}^^^!, 4>2\ - Rk)Pk[(t>2\ are obtained as the 
zeroth-orcler terms of the expansion of the modified prop- 
agator Pk-( a ,x )(6,a 3 ) that generalizes Eqs. (7) and (10). 
The above ERGE coincides with that previously derived 
without the superfield formalism by means of an expan- 
sion in number of free replica sums (when evaluated at 
T = 0). The same is true for the ERGE for all higher- 
order cumulants: for explicit expressions, see @. 

The (super)symmetries of the modified action in Eq. 
(5) are linearly realized and induce a set of WT identities 
for the 1PI generating functional T^,. 8 Taking functional 
derivatives of these identities with respect to the super- 
field and evaluating the resulting relations for superfield 
configurations $(x) = <fi(x) leads to relations for the cu- 
mulants. The most powerful relations mix cumulants of 
orders p and (p+1), the first nontrivial illustration of the 
latter kind being 



^rg> i;Ba [0 >( 
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which for fields that are also uniform in the Euclidean 
space gives a relation similar to that for the cutoff func- 

tions: rjgV;0,$ = A s^rg(g 2 ; <j>). 

An important feature of the present superfield theory 
is that SUSY leads to DR: this is obtained nonperturba- 
tively by combining the WT identity in Eq. (12) with the 
ERGE for the first cumulant in Eq. (11) and by follow- 
ing the line of reasoning of Refs. Q. As one knows that 
DR does not hold in low enough dimension, what then 
goes wrong in the formalism ? The answer is that SUSY, 
more precisely invariance under the superrotations when 
the theory is restricted to a single copy, is spontaneously 
broken along the flow and that a singularity occurs. From 
an analysis of the structure of the flow equations, we ex- 
pect that breaking of DR requires the presence of a linear 
"cusp" in the field dependence of rj^ , cusp that should 
appear at a finite scale during the RG flow. (On the 

other hand, weaker nonanalyticities in and nonan- 

alyticities in higher-order cumulants can only appear at 
the fixed point, in the limit k — > 0, thereby preserving the 
DR property.) This of course must be checked in actual 
calculations, which is what we provide below. 
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FIG. 1. NP-FRG flow of the dimensionless cumulant S k (f + 
y, <p — y) in d — 4 < don for <p = and for initial conditions 
close to the critical point. A linear cusp in \y\ appears at a 
finite RG "time" \t\ = log(A/fc). 



If SUSY is spontaneously broken, how can one continue 
the RG flow for the effective average action ? The original 
formal construction a priori loses its meaning, but a non- 
trivial continuation can be found if (i) one assumes that, 
except for the superrotations, all of the properties and 
symmetries of Tfe remain valid; most importantly, this 
includes the "Grassmannian ultralocality" encompassed 
in Eq. (10) that enforces single-solution dominance^, (ii) 
one only considers ERGE's for cumulants evaluated for 
generic (nonequal) field arguments, so that a putative 
nonanalytic dependence can freely emerge, and (iii) one 
modifies the regulator by replacing by a running 
which is the typical strength of the renormalized random 
field at scale k. More specifically, and in order to reach a 
fixed point and a scale- free solution describing the critical 
behavior of the RFIM, we choose R k (q 2 ) = Z k k 2 r(q 2 /k 2 ) 
and R k {q 2 ) = — (A k /Z k )r'(q 2 /k 2 ); Z k and A k are respec- 
tively obtained from d q 2T^{q 2 ) and T^^g 2 ) evaluated 
for q 2 = at zero field, and choices for the function r are 
given m ||,[l3 andOj. From Eq. (12) and below, one 
can see that so long as SUSY is not broken, Afc = AsZ k 
and the regulator is SUSY invariant, which guarantees 
the consistency of the RG description. 

Finally, we provide a SUSY-compatible nonperturba- 
tive approximation scheme for the ERGE. We combine 
truncations in the derivative expansion, which approxi- 
mate the long-distance behavior of the 1PI vertices, and 
in the expansion in cumulants of the renormalized disor- 
der. The WT identities require that the orders of trun- 
cation in the two types of expansions be related. The 
minimal truncation that can already describe the long- 
distance physics of the RFIM and does not explicitly 
break SUSY is the following: 



r fcl iM = 
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FIG. 2. Anomalous dimensions rj and 77 versus d. DR is ob- 
served above cLdr — 5.1. r? and 77 satisfy the required upper 
(77 < 2r;) and lower bounds (red dashed lines).— Stars corre- 
spond to simulation results . 15 ' 16 The region just below <1dr 
is unfortunately numerically difficult to access. 



with the higher-order cumulants set to zero. Inserted 
in the ERGE for the cumulants, the above ansatz pro- 
vides 3 coupled flow equations for the 1-copy poten- 
tial Uk{<j>) that describes the thermodynamics of the 
system, the field renormalization function Zk (</)), and 
the 2-copy potential Vft(0i,02) from which one obtains 
the second cumulant of the renormalized random field, 

rf^V = 0;4> u <h) = &k(<t>i,<h) = dfrdfrVkifafa). 

To search for the fixed point that controls the critical 
behavior (associated with a spontaneous breaking of the 
Z2 symmetry), the flow equations must be recast in a 
scaled form. The fixed point being a zero-temperature 
oneji^ the spatial decay of the correlations [see be- 
low Eq. (11)] at criticality is now characterized by two 
"anomalous dimensions" rj and fj, 

P{r) ~ r- (d - 2+,,) , P(r) ~ r^ 44 ^, (14) 



with 77 < fj < 2r), and one has to introduce scaling di- 
mensions involving an additional critical exponent. 6 The 
resulting equations are generalizations of those shown in 
Ref. and are not displayed here. We have solved these 
coupled partial differential equations numerically, look- 
ing for the proper (critical) fixed point as a function of 
dimension (more details will be given elsewhere). This 
procedure is numerically very demanding and requires 
handling 3 coupled equations for 2 functions of 1 vari- 
able (Uk and Zk) and 1 function of 2 variables (A&). 

An important property of the present theory is that if 
in the limit fa -> 0i, &k(4>i><h) = Afc O (0) + A fc2 (0)(0i - 
2 )2 _)_... w jth = 4- fa)/2, then the flow of A fe0 (^) 
coincides with that of Zk(4>): this is precisely the WT re- 
lation derived from Eq. (13), and DR exactly follows. On 
the other hand, a spontaneous breaking of the SUSY and 
of the associated WT identity occurs whenever Ak2(4>) 
diverges and A& has a cusp-like singularity in the form 

A*; (</>i, <?(>2) = Afco (</>) + Aka (</>)l0i _ <fo H as 02 -> 01- 

We find that the solution without cusp is stable and 
that 77(d) = 77(d) = r] Ism3 (d — 2), in agreement with the 
DR prediction, above a critical dimension dr)R ~ 5.1. 
For d < dun, we obtain a once unstable "cuspy" fixed 
point (see Fig. 1) and DR is broken: the exponents 77 and 
fj bifurcate, with 77(d) < fj(d) (see Fig. 2). In d = 3, we 
find rj ~ 0.57,7? ~ 1.08 and in d = 4, 77 ~ 0.24, fj ~ 0.40: 
this is in good agreement with the existing estimates; 15 ' 16 
which gives support to the whole scenario (the results are 
also 1-loop exact near d = 6). In addition, the continuous 
variation of 77 and 77 with d and the existence of a critical 
dimension above which 77(d) = 77(d) contradicts the claim 
that the two exponents are always related by a fixed ratio 
77(d) = 2r](d)H 

In conclusion, the present study provides key pieces for 
a complete resolution of the long-standing puzzles associ- 
ated with the critical behavior of the RFIM. In doing so, 
we have developed tools that may prove useful in other 
contexts where the need to select a unique solution of a 
stochastic field equation arises, as in "glassy" systems, 
turbulence or nonabelian gauge field theories. 
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